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Galaxy bimodality versus stellar mass and environment 
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ABSTRACT 

We analyse a z < 0.1 galaxy sample from the Sloan Digital Sky Survey focusing on the 
variation of the galaxy colour bimodality with stellar mass Ai and projected neighbour 
density S, and on measurements of the galaxy stellar mass functions. The character- 
istic mass increases with environmental density from about 1O 1O 6 A^0 to 10 10 ' 9 .M© 
(Kroupa IMF, H = 70) for S in the range 0.1-10 Mpc -2 . The galaxy population 
naturally divides into a red and blue sequence with the locus of the sequences in 
colour-mass and colour-concentration index not varying strongly with environment. 
The fraction of galaxies on the red sequence is determined in bins of 0.2 in log E and 
logVW (12 x 13 bins). The red fraction f r generally increases continuously in both E 
and Ai such that there is a unified relation: f r — F(E,Ai). Two simple functions are 
proposed which provide good fits to the data. These data are compared with analogous 
quantities in semi-a naly tical models based o n the Millennium TV-body simulation: the 
iBower et all (120061) and lCroton et all {2006) models that incorporate AGN feedback. 
Both models predict a strong dependence of the red fraction on stellar mass and en- 
vironment that is qualitati vely similar to the observations. However, a quantitative 
comparison shows that the IBower et alJ model is a significantly better match; this 
appears to be due to the different treatment of feedback in central galaxies. 



Key words: galaxies: evolution 
luminosity function, mass function. 
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1 INTRODUCTION 

Galaxies when characterised by their morphology or ra- 
dial profiles, integrated or central colours, and total lu- 
minosity or stellar mass, exhibit a range of rel ationships. 
These include colour-mor phology relations iHolmberelll958t 
iRoberts fc Havnesl Il99l) . and colour-magnitude relations 
separatel y for earlv-tvpe galaxies <TFaberll973l) and late-type 
galaxies dChester fc Robertsi ri964'l. While it was often con- 
sidered that the natural dividin g line was between spirals 
and ellipticals/lenticulars ("e.g. iTullv. Mould fc Aaronsonl 
1982), it was not until the multi- wavelength Sloan Digital 
Sky Survey (SDSS) that the galaxy population was con - 
sidered strongly bimodal in colour dStrateva et alJ l200ll) . 
Even when considering other galaxy properties such as ra- 
dial pr ofiles, the natural division is into two galaxy popu - 
lations llHogg et alJl2002t lEllis et~al]l2005t iBall et alJl200r3h 
Given that large automated imaging surveys are better at 
defining a galaxy's colour than morphology, it is more nat- 
ural to describe a galaxy as being on the "red sequence" 



or "blue sequence" rather than being an "early type" or 
"late type" . This interpretation also has the advantage that 
galaxy colours are directly related to the star formation, 
dust and metal enrichment history of the galaxy and can 
thus be more readily interpreted in theoretical models. 

A key goal of galaxy evolution theory is to ex- 
plain the bimodality and the relationships within each se- 



recently (e.g. |Kan£_et_alJ |2005j; | 


Menci et alJ l200l 


120061 


Springel et al J 2005b 


Avila- Reese 


2006; 


Bower ct al. 


2006; 


Cattanco ct al. 2006; 


ICroton et al 


2006; 


Dckcl & Birnboirn 


2006; Perez ct al. 2006J. The key ingredient of many of these 



models is the inclusion of feedback from active galactic nuclei 
(AGN). Although AGN feedback (or its equivalent) is imple- 
mented in different ways in each of these models, the overall 
e ffect is to suppress c ooling in massive halos. For example, 
in lBower et alJ i2006T) the AG N feedback suppress es 
hydrostatic cooling flows, while ICroton et alJ (I2006T) adopt a 
semi-empirical description for the AGN power related to the 



2006) the AGN feedback suppresses quasi- 
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Bondi accretion rate. Bimodality of galaxy colours does not 
directly result from these schemes, however, the AGN feed- 
back allows the star formation rate parameterisation to be 
adjusted in such a way as to simultaneously obtain a good 
description of the colour distribution at faint magnitudes 
and a good match to the shape of the luminosity function. 
In these models red galaxies at faint magnitudes are pre- 
dominantly satellite galaxies of brighter systems, while at 
bright magnitudes the central gal axies are also red because 
of the AGN feedback (e.g. fig. 3 of lBower et all) . In the real 
universe, however, the association between galaxy colours 
and their location in t he halo is unlikely to be so simple 
JWeinmann et aljEoOrj . and measurements of the depen- 
dence of galaxy colours on luminosity and redshift are an 
important constraint on the new generation of galaxy for- 
mation models. Analysis of the relationship between galaxy 
colours and environment will place important constraints on 
the processes defining galaxy evolution. 

This paper provides a detailed analysis of the varia- 
tion of the bimodality with stellar mass and environment. 
The plan of the paper is as follows: in § 11.11 we review 
a series of papers providing the buildup to this paper; 
in §|2] we describe the data; in §[3] we present the re- 
sults; in §2] we discuss the implications and compare the 
data with models; and in §|5| we summarise the main re- 
sults. Fits to the mass functions are presented in the Ap- 
pendix. The data represented in this paper is available at 
http://www. astro. livjm.ac.uk/^ikb/research/ or upon re- 
quest. 

1.1 Previous work 

iBaldrv et alJ i2004al hereafter Paper I) characterised the 
volume-averaged colour-magnitude distribution of galaxies 
by fitting double Gaussian functions to the colour his- 
tograms in magnitude bins. This showed that the red se- 
quence transitioned from a broader colour distribution at 
M r ~ —19 to a narrow distribution for massive early 
types at M r ~ —21 while the blue sequence become 
significantly redder over the range —20 to —22 (using 
centrally- weighted colours). This supported the suggestion 
of a transition in galax y properties around 10 10 ' 5 Mq 
jKauffmann et alJ [2003b) and the former is consistent 
with faint red-sequence galaxies forming more recently 
iDe Lucia et al]|2004t lKodama_et_aL| ]2004|j_but not neces- 
sarily in the richest clusters lfATcteonTl200^) . The differ- 
ence in star formation history for galaxies below and above 
the transition mass may be relate d to the balance between 
hydrostatic versus rapid coolin g teirnboim^^JOekell [2 0031: 
iKeres et al1l2005l: iBower et aljEobet ICroton et al.ll2006t) ."or 
to the diversity of star formatio n histories of satellite galax- 
ies at low masses (lBowerlll99ih . 

Balog h et alJ (I2004al . hereafter Paper II) analysed Ha 
emission strength as a function of galaxy environment. The 
Ha equivalent width (EW) distribution is bimodal with 
the distribution of the star forming population (blue se- 
quence) not depending strongly on environment. The frac- 
tion of galaxies with EW > 4A varied continuously with 
environmental density and there was no evidence of a break 
density that had b een reported before llLewis et al l b mi 
iGomez et al.l I2003T) . This demonstrated the importance of 
describing the variation of the bimodal population by com- 



paring the number of galaxies within each population rather 
than using a quantity averaged over both populations. 

In Paper II, the results indicated a dependence of the 
star- forming fraction on scales of about 5Mpc (after ac- 
counting for local environment). However, the interpretation 
is difficult because the small-scale measurement is noisy and 
the large-scale measurement could actually b e adding infor- 
matio n ab out the small-sca l e env ironment felanton et alJ 
2006), and lKauffmann et alJ (120041) found no environmental 
dependence on scales larger than 1 Mpc. It is also known 
that for galaxies in clusters properties can depend on en- 
vironment for scales smal ler than 1 Mpc (iDressled Il98d : 
IWhitmore fc Gilmoreil99lh . These argue for an iVth nearest 
neighbour approach to measuring environmental density as 
the radius is smaller in high density regions and expands in 
low density regions where there are often insufficient galaxies 
on s mall scales (for a fixe d luminosity cut). 

iBaloeh et alJ J2004bl . hereafter Paper III) extended the 
double Gaussian fitting of Paper I to colour histograms 
across environment and luminosity. At fixed luminosity, the 
mean positions of the sequences become marginally redder 
with environmental density. For the red sequence, this can 
be explained by a small difference in a ge between low- and 
high-de nsity environments of ~ 2 Gvr ([T homas et alJl20o3) 
or less (iHogg et al]l2004t liBernardi et alJl2006l) . In contrast, 
the fraction of red-sequence galaxies varied strongly with 
environment as measured by projected density. No effect re- 
lated to velocity dispersion of a group or cluster was detected 
which is consistent with the morpholo gy-density relation be - 
ing similar in groups and clusters llPostman fc Gelleilll984h . 

IBaldrv et al] ]2004bl) confirmed the shifts in colour of 
the red and blue sequences, 0.05 and 0.1 in u—r, respectively, 
over a factor of 100 in projected density. The red fraction 
varied from 0% to 70% for low- luminosity gal axies and 50% 
to 90 % for high- luminosity galaxies (see also iTanaka et al] 
2004). The effects of environment and luminosity could be 
unified in that the fraction of red-sequence galaxies was 
related to a combined quantity: S mo d = (E/Mpc~ 2 ) + 
(L r /L rjnorm ) where E is the projected density and L rjnorm 
is luminosity of a galaxy with M r — —20.2. In this paper, 
this effect is explored in more detail using a larger sample 
and by converting luminosity to stellar mass. 



2 DATA 

In this section, the selection of data from the SDSS catalogue 
and the derived quantities for each galaxy are described. The 
basics of the SDSS are described in § l2.1l while the sample se- 
lection is described in § 12.21 The primary sample consists of 
spectroscopically observed galaxies but a larger sample was 
used to determine spectroscopic completeness (§ 12.311 and to 
measure uncertainties in environmental densities (§ 12.611 for 
which photometric redshifts were determined (§ 12.41 . Num- 
ber density corrections are described in § 12.71 The derived 
quantities include absolute magnitude M r , rest-frame colour 
Cur, inverse concentration index C (§ 12.511 . projected density 
of neighbouring galaxies E (§ 12.611 . and an approximate stel- 
lar mass M (§ 12.81 . 
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Figure 1. Sky coverage for SDSS DR4. The solid lines outline the primary photometric area, the dashed lines outline the galactic poles at 
20° latitude, and the points show the positions for 4% of the galaxies that have been observed spectroscopically. The edges are are easier 
to define in SDSS survey coordinates, which is a spherical coordinate system [(eta,l ambda)=(0,0) corresponds to (r a,dec)=(185,32.5)]. 
The photometric area covers 6670 deg 2 and the spectroscopic area covers 4780 deg 2 (Add man- McCarthy et al.H2 006 ) . 



2.1 Sloan Digital Sky Survey 

The SDSS is a project, with a dedicated 2.5-m telescope, 
desig ned to image 10 4 deg 2 and obtain spec tra of 10 6 ob- 
jects jYork et allboOOHStoughton et alJl2002ft . The imaging 
covers five broadbands, ugriz with effective wavelengths of 
355, 467, 616, 747 and 892 nm, using a mosaic CCD camera 
ijOunn et alJIlQfiisF) . Spectra are obtained using a 640- fibre 
fed spectrograph with a wavelength range of 380 to 920 nm 
and a resolution of A/AA ~ 1800. 

The images are reduced using a pipeline photo that 
measures the observing conditions, and detects and mea- 
sures objects. In particular, photo produces various types 
of magnitude measurement including: (i) 'Petrosian', the 
summed flux in an aperture that depends on the surface- 
brightness profile of the object; (ii) 'model', a fit to the flux 
using the best fit of a de-Vaucouleurs and an exponential 
profile; (iii) 'PSF', a fit using the local point-spread func- 
tion. Th e magnitudes are extinction-correcte d using the dust 
maps of lSchleeel. Finkbeiner <S ^Davij (|r^98|). Details of the 
imaging pipelines are given bv lStouehton et al.l (12002ft . 

Once a sufficiently large area of sky has been imaged, 
the data are analysed using 'targeting' software routines that 
determine the objects to be observed spectroscopically. The 
main sets of targets are: (i) the main galaxy sample (MGS), 
an r < 17.77 selection iStrauss et al J 120021: (ii) luminous 
red gal axies (LRGs), an r < 1 9.5 and colour-selected galaxy 
sample ( Eisenstcin ct al 1200 lft : and (iii) quasi-stellar objects 
(QSOs: lRichards et al.ll2002ft . The targets from all the sam- 
ples are then assig ned to plates, each wi th 640 fibres, using 
a tiling algorithm jBlanton et al.ll2003aft . The main restric- 
tion is that two fibres cannot be placed within 55" on the 
same plate. 

Spectra are taken using, typically, three 15-minute ex- 
posures in moderate conditions (the best conditions are used 
for imaging). The signal-to-noise ratio (S/N) is typically 10 
per pixel (~ 1-2A) for galaxies in the MGS. The pipeline 
spec2d extracts, and flux and wavelength calibrates, the 
spectra. The spectra are then analysed by another pipeline 



that classifies and determines the redshift of the object. The 
redshift success is about 99% for galaxies considered in this 
paper. 

2.2 Sample selection 

The data use d in this paper are from the SDSS Data Release 
Four (DR4) iAdelman-McCarthv et"aill2006ft . The primary 
sample consists of 151 642 galaxies with r < 17.77 in the red- 
shift range 0.010-0.085 (Sample C), with environmental den- 
sities determined from the distance to the 4th and 5th near- 
est neighbours with M r < —20. However, the spectroscopic 
sample is not 100% complete; some galaxies are missed be- 
cause of fiber-placement restrictions or are not yet observed, 
and there are edges to the sky coverage. Figure^shows the 
sky coverage for the photometric and spectroscopic data. In 
order to accurately account for missed spectra and to calcu- 
late completeness corrections, a larger photometric sample 
was considered. The galaxy sample selection and calculation 
of environmental densities are outlined below. 

Objects were selected from the 'CasJobs' catalogue 
archive server 1 with criteria 10.0 < r < 18.0 , where r is 
the Petrosian magnitude corrected for Milky-Way extinc- 
tion, and rpsF — r mo dei > 0.25 . Photometric data ('Galaxy' 
table) and spectroscopic data ('SpecObjAU' table) where 
available were extracted. This produced an initial sample of 
1016 565 objects after removing duplicates (Sample A). 

The next stage was to remove low-surface brightness ar- 
tifacts and stars. A modified surface brightness was defined: 
/timod = Mr,50 — 0.3(r — 18) where [i r ,50 is the mean sur- 
face brig htness (SB) within the Petrosian half-light radius 
(eq. 5 of IStrauss etai]l2002ft . The slope of 0.3 per magni- 
tude was determined empirically from the SB-magnitude 
distribution so that the modified quantity was a better 
separator for stars, galaxies and artifacts. Figure [5] shows 
SB versus magnitude for various samples and the dividing 

1 http:/ /casjobs. sdss.org/CasJobs/ 
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Figure 2. Surface brightness versus magnitude. The grey points 
and dashed contours represent all objects in Sample A; the red 
points and contours represent spectroscopically confirmed galax- 
ies; and the blue points and contours represent stellar systems. 
The green lines show the limits imposed for photometric objects 
selected for Sample B. Note that most stars are rejected by the 
'not saturated' criteria rather than the lower limit shown here. 



lines. Objects with no spectroscopy were then trimmed to 
18.5 < Hmod < 24.0 and 'not saturated' while spectroscopic 
objects were selected with 0.005 < z < 0.3. The trimmed 
sample contained 850 906 galaxies (Sample B). 



2.3 Spectroscopic completeness 

Spectroscopic completeness was determined from a related 
sample that included all spectra and was reduced to spec- 
troscopic regions (634 395 galaxies, 478 385 with r < 17.77). 
Galaxies were divided into classes based on their magnitude 
and the number of r < 17.77 neighbours within 55" (0, 1 or 
> 2). The completeness was given by the number with spec- 
tra divided by the total number in a class. Figure |3] shows 
completeness versus magnitude. As expected the complete- 
ness is low above 17.77 and it becomes colour dependent 
(LRG selection). Thus, the final sample (Sample C) is re- 
stricted to this MGS limit where the completeness correc- 
tions are reliable. It is also clear that in order to determine 
the relative numbers of galaxies in different environments, 
it is important to account for the change in completeness as 
a function of the number of neighbours. 



2.4 Photometric redshifts 

Photometric redshifts were determined for Sample B using 
X 2 fitting to 'training set' galaxies with spectroscopic red- 
shifts. The training set consisted of 396 302 galaxies out of a 
possible 410 206, with some galaxies rejected because of low 
redshift confidence, or extreme colours or colour errors. The 
data fitted for each galaxy were rp e tro and 4 'model' colours: 
u ~ 9> 9 ~ r i r ~ h i ~ z - Thus, the underlying assumption 
is that the redshift-magnitude-colour distribution is simi- 
lar for the galaxies with missing spectra as for the training 
set. Softening errors were added in quadrature, prior to de- 
termination of the % 2 values. These were 0.1 for the mag- 



Figure 3. Spectroscopic completeness versus magnitude. The 
three lines represent galaxies with 0, 1 or 2 neighbours within 
55" (corresponding to the 'fibre collision' radius). The arrow 
shows the limit for MGS selection; the completeness drops prior 
to that because of changes in catalogue photometry and targeting 
algorithms. Here we have used the DR4 best photometry. 

nitude (larger for r < 13) and 0.02 for the colours. 2 The 
photometric redshift is then given by the weighted mean of 
training-set redshifts that fitted within x 2 < 15 (99% limit 
for 5 degrees of freedom) . The weight for each match is given 
by exp(— Xi/2) w(zi), where w{z%) is a minor adjustment to 
reduce the implicit weight given to redshift peaks within 
the training set. Objects with no matches within the limit 
or with few matches in the range 10 < x 2 < 15 were as- 
sumed to be outside the redshift range of the training set 
(i.e. stars, or z > 0.3) and not considered in further anal- 
ysis (13 741 objects). An initial photometric redshift error 
was determined from the standard deviation (with weights) . 
These errors were then multiplied by a factor to give 95% 
confidence-limit errors: 

r ( l^photo ^spcc| \ . , 

Zerr,95 — /95 I I ^crr, initial 1,1) 

\ -2 crr) initial / 

where /gs() is the 95-th percentile of the expression deter- 
mined for galaxies with spectroscopic redshifts. After cali- 
bration, the median value of Zerr,95 is 0.04 with typical val- 
ues from 0.02 to 0.07. A minimum error of 0.02 was used in 
subsequent analysis. 

2.5 Absolute magnitude, rest-frame colour and 
concentration index 

The r-band absolute magnitude used in this paper is given 
by 

M r = r-fc r -51og(I3i/10pc) (2) 

where: r is the Milky-Way-extinction-corrected Petrosian 
magnitude; Dl is the luminosity distance for a cosmology 

2 The x 2 values were determined only using the errors for the 
galaxy being considered, not the training-set galaxies, otherwise 
higher weight would be given to training-set galaxies with larger 
errors! For testing and calibration of the redshift errors, training- 
set galaxies were not matched to themselves. 
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with (fi m ,fi A )o = (0.3,0.7) and H Q = 70 kms-'Mpc" 1 
and; k r is the k-correction using the method of lBlanton et alJ 
( 20023) ; 3 and rest-frame colour considered is given by 



C ur (t^modcl ) ('"model ^-r) 



(3) 



This was chosen as the optimum bimodality divider (Pa- 
per I). In effect, 'model' colours are centrally- weighted 
colours; the profile used is the best fit defined in the r-band 
light. 

Galaxies with photometric redshifts have a range of un- 
certainty in their absolute magnitude. For calculation of 
environmental densities, their fc-corrections were approxi- 
mated by k r = [0.3 + 0.3(u - r) mo del]ztest where zteet takes 
different values: the spectroscopic redshifts of the galaxies 
for which environmental density is measured. 

The (inverse) concentration index was used in some of 
our analysis. The is given by C — R50/R90 where R50 and 
R90 are the radii containing 50% and 90% of the Petrosian 
flux, averaged for the r and i bands. For typical galaxies, 
C ranges from 0.3 (concentrated) to 0.55; for comparison, a 
uniform disk would have C = 0.75. 



2.6 Environmental densities 

Densities were determined for Sample C using the infor- 
mation from Sample B. For environmental densities around 
each galaxy, we determine 



Ew = 



N 
nd 2 N 



(4) 



where cLm is the projected comoving distance to the Nth. 
nearest neighbour that is a member of the density defining 
population (DDP) and that is within the allowed redshift 
range. The redshift range is ±Azc = 1000 km/s for neigh- 
bouring galaxies with spectroscopic redshifts and ±z er r,95 
for those with photometric redshifts only. The DDP are 
galaxies with M r < M r ji m it — Q(z — zo) w here Mrjimit = 
-20, Q = 1.6 is the evolution determined bv lBlanton et all 



(2003c) and zo = 0.05. For galaxies with photometric red- 
shifts only, M r was determined as if the galaxy were at the 
redshift of the galaxy whose environment was being mea- 
sured. If the distance to the photometric edge d c d go were 
less than djv then nd 2 N in Eq. 21 was replaced by the area 
within a chord crossing the circle, i.e., assuming the edge 
was straight. 

The best estimate E was obtained by averaging log E m 
for TV = 4 and 5, and by averaging the calculation with 
spectro-z only and the full sample B. In addition E 

mill WaS 

determined using djv with spectro-z only, and E max using 
the smaller of djv with Sample B and d c d gc . The uncertainty 
in log E was given by the larger of log E max — log E and 
log E— log Emiii. For various environmental analyses, galaxies 
were only included if the uncertainty were less than 0.4 dex 
or less than 1 bin for galaxies in the minimum or maximum 
density bins. This more complicated procedure than simply 
using spectro-z only and discarding galaxies near an edge 
ensures that the maximum number of galaxies can be used in 
any analysis; and it accurately accounts for 'fibre collisions'. 



3 The fc-corrections were derived from kcorrect v. 4.1.3. 
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Figure 4. The number densities of galaxies versus redshift. The 
squares with error bars (including photometric redshifts) and di- 
amonds (spectroscopic redshifts only) represent the number den- 
sities of galaxies more l uminous than —20 (ev olution corrected to 
z = 0.05 with Q = 1.6. lBlanton et al.lEu03cl) : this is the density 
defining population used in this paper. The crosses with error bars 
and triangles represent galaxies more luminous than —21. Above 
about z > 0.09 and z > 0.13, the counts become incomplete in 
the r < 18 sample for the M r < —20 and M r < —21 popula- 
tions, respectively. At z < 0.02 there is incompleteness because 
of problems with automatic determination of photometric prop- 
erties. Using Q = 1.6 for the luminosity evolution, the number 
density evolution is near flat across the reliable redshift range, 
shown by the dashed lines. 



Variations in the above parameters (Azc = 1000 km/s, 
A/rjimit = —20, N = 4 and 5) used to determine E were 
tested. The results were not strongly dependent on the exact 
values, and the values chosen were near optimum and consis- 
tent with earlier papers (Paper II and Paper III). For typi- 
cal galaxies, E ranges from 0.05 Mpc -2 (void-like regions) to 
20 Mpc -2 (clusters); for comparison, E = 0.5 Mpc -2 for the 
median environment around DDP galaxies whereas E would 
be 0.14 Mpc" 2 if DDP galaxies were distributed uniformly. 

Figure 01 shows the number density of the DDP popu- 
lation versus redshift and the density of a more luminous 
Mr < —21 population. The near constant number density 
up to the completeness li mit justifies the use of th e evolution 
parameter determined bv lBlanton et all i2003c|) . Sample C 
was restricted to the redshift range 0.010 to 0.085 for relia- 
bility of the environmental density measurements. 



2.7 Volume and density corrections 

The DDP is volume limited but we also consider galaxies 
fainter than M r — —20. For these and for some galaxies near 
—20, a volume correcti on is applied to the number density 
using a V max method (lFeltenlll976l^ . Here, the weight of a 
galaxy for computing number densities is given by a factor 
1/Vmax where Vmax is the maximum volume over which the 
galaxy could be observed within the r < 17.77 limit and 
the survey volume (0.010 < z < 0.085 over 1.52 steradians; 
2.3 x 10 7 Mpc 3 ). 

The added complication when considering environment 
is the fact that the contribution of a particular environment 
bin to the volume-averaged number density can vary sig- 
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nificantly as a function of redshift. To correct for this, we 
determine a relative volume (fv) that a galaxy in a partic- 
ular environment could be observed over the redshift range 
given by Vm ax (0.010 < z < z max ). This relative volume is a 
function of z max and E and can be determined by measuring 
the number densities of bright galaxies that can be viewed 
over the whole survey volume. The relative volume was de- 
termined in bins of 0.05 in z and in E (depending on the 
bins being used). In other words, fv is the number density 
of bright galaxies with z < z max relative to the number den- 
sity covering the whole redshift range (with reduced survey 
area to avoid edge effects). The number density of galaxies 
with a given 2 max and E was renormalised by 1/fv- This 



factor fv is in the range 0.6-1.1 for 



> 0.03 over 12 



environment bins; fv is on average less than unity because 
galaxies are rejected when the uncertainty in E is too large. 

To summarise, there are three number density correc- 
tions for each galaxy in Sample C: (i) 1/C S where C 3 is the 
spectroscopic completeness which is a function of r and the 
number of MGS neighbouring galaxies within 55" (S I2.3H : (ii) 
1/Knax where V max is a function of the absolute magnitude 
and the galaxy's colours which determine fc-corrections; (iii) 
1/fv where fv is the relative volume which is a function of 
z max and E when dividing the population into environment 
bins. 
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Figure 5. Correlation of galaxy stellar mass-to-light ratio with 
colour. Stellar masses were determined for a sample of 10 4 galax- 
ies. The blue squares and dotted-line region represent the me- 
dian in 0.1 color bins and 1-sigma ranges computed by K.G. 
jGlazebrook et alj EooA. The green diamonds and region rep- 
resent stellar masses co mputed by the MPA Garching group 
iKauffmann et aljEooSah . The red dashed line represents a bro- 
ken line fit to the data, with the coefficients shown in the plot. 



2.8 Conversion to stellar mass 



3 RESULTS 



widely promoted (e.g. iTh 


ronson & Greenhouse! 


198S; 


Prueniel & Simienl Il99rl 


Brinchmann & Ellis 


2000; 


Benson et al.l 120021: iBell et al 


120031 iPerez-Gonzalez et alJ 



Parmelia^^nTl2006l) . The total stellar mass of a galaxy 



is a more fundamental quanti ty than luminosit y as ha s 
been shown by, for example, IKauffmann et al] (l2003aT) . 
Galaxy stellar masses are typically derived by fitting 
stellar-population synthesis models to colours or spectra. 
However, there are systematic uncertainties that depend 
on the underlying set of allowed star-formation histories in 
the models as well as uncertainties in th e models. Here, 
we us e an approach similar to that of IBell fc de Jond 
(l200lt) whereby we use a stellar mass-to-light ratio that 
is a function of one colour only. While uncertainties will 
remain, the method has the advantage of retaining a simple 
relation between the derived physical quantity, the selection 
function, and the rest-frame colour used in our analysis. 
Figure |5] shows a plot of stellar mass-to-light ratio versus 
Cur colour for 10 4 galaxies analysed by the MPA Garching 
group and by K.G. The dashed line represents the relation 
that we use in our analysis. 4 

4 We use M for mass and M for absolute magnitudes. The con- 
version is given by \o r(M/Mp)) = (M rW — M r ) /2.5 + log(A4 / L r ) 
where M rQ = 4.62 <Blanton et al.ll200lT) . The computation of 
stellar masses by the MPA Garching group uses the spectral fea- 
tures, D4000 a nd H<5, and the Pe t rosian z-band magnitude and 
is described by IKauffmann et al] <2003a|): while the method of 
K.G. is described by iGlazebrook et alJ fedolh except Petrosian 
magnitudes were u sed here . For b oth analyses, JVl/L r ratios were 
obtained using the lKroupal fcOOll) IMF (eq. 2 of that paper, from 
O.IMq to 100 or 120At Q ). 



3.1 Colour-concentration relations as a function 
of stellar mass 

In Paper I and Paper III, we divided the galaxy population 
by using double Gaussian fitting to the colour functions in 
luminosity bins. This separates the galaxy population into 
red and blue sequences. There was no consideration of the 
morphology or structure of the galaxies. This method is ac- 
curate if the colour functions are truly the sum of Gaus- 
sian distributions and if the errors are small. A more robust 
method of outlining the two sequen ces is to consider a joint 
distribution in colour and structure JPriver et al.l200rj) . Fig- 
ure H2 shows the distribution of observed galaxies in colour 
versus concentration index for three different stellar mass 
ranges. 

The peak of the red sequence gets redder (2.3 to 2.6) 
and more concentrated (0.38 to 0.32) with increasing mass 
(logA-1 from 9.5 to 11.0); while the blue sequence gets red- 
der (1.5 to 2.1) with approximately the same concentration 
(0.45). This means that the natural dividing lines of the pop- 
ulation are around C = 0.4 for all the stellar masses (vertical 
dashed lines in the figure) but the dividing line in colour gets 
redder (horizontal dashed lines). At lower masses, dividing 
using colour is significantly better than dividing using con- 
centration index. At higher masses, while colour is still a 
reasonable divider, concentration index can be useful in ro- 
bustly determining the colour mean and dispersion of each 
sequence. 

In Paper I, a method for fitting double Gaussians to the 
colour functions is described. At each luminosity or mass 
bin, there are six parameters, amplitude, mean and disper- 
sion for each sequence, given by <p T , /Lt r , 0> , 4>b, A*b, When 
the counts are low or the sequences are significantly merged 
together in colour (e.g. at the high-mass end), the solution 
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Figure 6. Colour versus concentration index for different stellar-mass ranges. The points and solid contours represent galaxies in our 
sample. The contours are logarithmically spaced in number density with 4 contours per factor of 10. The dashed lines show the position 
of possible colour or concentration-index dividers. 



is not well defined. With the robust method, /i r and a v are 
determined using galaxies with predominately high concen- 
tration (full weight for C < 0.4) while ^ and <Tb are deter- 
mined using galaxies with low concentration (C > 0.4). The 
amplitudes of each Gaussian are still determined using a fit 
to both populations but with the /j, and a values fixed. 

This method of determining the positions of the se- 
quences is used in § 13.21 which then determines the best-fit 
dividing line in color, and in § 13.41 In § 13.51 only the dividing 
line is used. 



3.2 Volume-averaged colour-magnitude and 
colour-mass relations 

Figure EJshows a comparison between colour-magnitude and 
colour-mass relations: C ur versus (a) it-band absolute mag- 
nitude, (b) r-band absolute magnitude, and (c) logarithm of 
the stellar mass. Note how the contours are almost vertical 
at the luminous end in Fig. EI a) whereas the red sequence 
is dominant at the high mass end in Fig. [TJ C ) ■ 

In order to fit the mean positions of the sequences, we 
used a tanh plus a straight-line function as per Paper I. A 
general form of this 'T function' is given by 



T(x) — po + pix + qo tanh 



'12 



(5) 



where po,i are the straight line parameters; and go, 1,2 are the 
tanh parameters. The fitting of the sequences followed the 
iteration outlined in § 4.2 of Paper I except with the addition 
of the robust method outlined in § 13. II above for the double 
Gaussian fitting. 

The grey regions in Fig. show the mean position and 
dispersion along each sequence. A key point is that the blue 
sequence is significantly narrower versus stellar mass than 
luminosity, particularly near-UV luminosity (the color dis- 
persion is about 0.23, 0.31, 0.37 for blue-sequence galaxies 
with log M ~ 10.5, Mr ~ -21.0, M u ~ -19.5, respectively). 
This demonstrates that stellar mass is closer to a fundamen- 
tal predictor of a galaxy's colour than luminosity. The dash- 



and-dotted lines in the figure represent the best-fit dividing 
lines between the red and blue sequences using the optimal 
criteria of Paper I. 

When dividing the galaxy population by environment, 
it becomes more difficult to proceed with the full fitting 
procedure. For cases where many environmental bins are 
used, the number of galaxies in each sequence is defined 
by galaxies redder and bluer than the best-fit dividing line 
[dash-and-dotted line in Fig. Etc)]: C ur ,divide varies from 1.8 
to 2.5 for log M from 9.0 to 11.5; T function with parameters 
ofpo.i = (2.18,0), 90,1,2 = (0.38,10.26,0.85). 



3.3 Environmental dependence of the galaxy 
stellar mass functions 

Before considering the variation of the colour-mass relations 
it is important to note that the galaxy stellar mass functions 
(GSMFs) vary significantly with environment. The galaxy 
population was divided into 12 environment bins and 16 
mass bins. Figure |H| shows the variation of the GSMFs with 
environment. The GSMFs are plotted using mass fraction 
normalised by total mass within each environment bin (i.e. 
the integral under each curve is unity). This has two advan- 
tages: (i) the t/-axes values range over 2 orders of magnitude 
compared to about 4 when plotting number density; and (ii) 
the peak of each GSMF corresponds to the dominant con- 
tribution to the total stellar mass within an environment 
bin. This varies from log M ~ 10.6 to 11.0 from low to high 
densities; with the fractional contribution due to high mass 
galaxies (log M ~ 11.5) increasing by a factor of about 50. 

The GSMFs can be characterised by fitting ISchechteil 

lll976ft functions: 



oc M a e~ M/M * 



(6) 



where 4>n &M is the number of galaxies with masses between 
M and M + dM ; or 

(pio g M oc Ai e (7) 
where 0i og m d log M is the total mass of galaxies between 
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Figure 7. Comparison between galaxy colour-magnitude and colour-mass relations. Colour is plotted versus (a) M u , (b) M r and (c) 
log Jv[. The dashed and solid lines represent logarithmically-spaced density contours with 4 contours per factor of 10; the number densities 
are completeness and volume corrected. The dotted lines represent galaxies with stellar masses of 10 10 and 10 11 Mq. The grey regions 
represent the colour means and ±l-sigma ranges of the red and blue sequences. The dash-and-dotted lines represent the best-fit dividers 
between the sequences based on the criteria of Paper I. 



1.00 



CO 



0.10 



0.01 




1.00 



0.10 



0.01 



9.5 10.0 10.5 11.0 11.5 
(M / -tt S oiar) (stellar mass) 




loi 



9.5 10.0 10.5 
(M / ^ S oiar) (stellar 



11.0 11.5 
mass) 



Figure 8. Galaxy stellar mass functions. In panels (a) and (b), the symbols and lines represent different environments as shown in the 
legend of panel (a). Note the unusual y-axis, with respect to standard mass or luminosity functions, shows clearly that the 'peak' in the 
mass function is near the Schechter break (log.M*). Panel (a) shows the binned functions (the median errors are 5%) while panel (b) 
shows the Schechter function fits. Inset panel (c) shows log A')* versus logS. The diamonds represent log Ai* as defined by the parameter 
in the fit while the squares represent the peaks of the fits given by log[A^*(« + 2)] . Note that a standard Schechter function fit is only 
valid over the range from about 9.5 to 11.5, and there is a statistically significant upturn in the faint-end slope at low masses (cf. Paper I, 
iBlanton et all2005bl Appendix) . Inset panel (d) shows the approximate fractional contribution to the stellar mass density of the universe 
(f2 sta rs ~ 0.002) for each environmental bin from high to low densities (top to bottom). Plots of the GSMFs using double Schechter fits 
and using number versus log M are shown in Fig. IA1I 
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Figure 9. Colour versus stellar mass relations for different environments. Panel (a) represents void-like environments while panel (f) 
represents cluster-like environments. The dashed and solid lines represent logarithmically-spaced density contours. The grey regions 
outline the red and blue sequences. The dash-and-dotted lines show the best- fit divider from Fig. ITTc - ). 



log.M and logAI + dlogA4. Normalised 4>\ os m is plotted 
in Fig. |H] For a standard Schechter function, the peak of 
the 4>i gM GSMF is given by X pC ak = M*(a + 2). Thus, 
A^p ea k can be considered an alternative characteristic mass 
for faint-end slop es around — 1, which is less dependent on 
model choice fcf. lAndreonll2004l) . 

The Schechter fits are shown in Fig. |SJb) with logAI* 
versus environment plotted in Fig. |HJc) . A fit to log -M pC ak 
is shown by the dotted line, formally 

log[AT(a + 2)] = (10.73 ±0.01) + (0.15 ± 0.01) log E (8) 

The faint-end slope is approximately —1 over the fitted range 
which does not vary strongly with environment. A standard 
Schechter function does not fit the entire mass range. Fits to 
the extended mass range using a double Schechter function 
(Paper I) are shown in the Appendix. 

The increase of the characteristic mass with local den- 
sity of 0.4 dex, equivalent to lmag, is broadly consistent 
with other results. A precise quantitative comparison is dif- 
ficult because of the different method s of determin i ng en - 
vironment and different wavelengths. ICroton et al.l <l2005h 
found an increase of ~ 1 mag betwee n void and cluster- 
like environments using bj luminosities. iHovle et alJ (l2005ll 
found a difference of 0.9 mag between void and wall galax- 
ies using r-band luminosity. There are many other analyses 



looking at luminosity fu nctions versus group masses, clus- 
ter properties, etc. (e.g . |Bingggli^^ilJ ^ggj [Balggh e t alJ 
l200lt iDe Propris et alJ l2003t IZandivarez et al.ll200rj) . The 

main point is that the overall characteristic luminosity or 
mass increases with environmental density. This should be 
accounted for when comparing galaxy colours between dif- 
ferent environments. 



3.4 Environmental dependence of the colour-mass 
and colour-concentration relations 

The galaxy population was divided into environmental bins 
in order to determine the variation of the colour-mass rela- 
tion. Figure |5] shows the colour-mass relations for six envi- 
ronmental bins ranging from log E < —0.8 to log E > 0.8. 
The equivalent of the cluster red sequence is evident in 
Fig. I^Jf) ; while the dominance of the blue sequence at lower 
masses is evident in Fig.^a-c) for the lower densities. While 
the red sequence is dominant at higher masses in all envi- 
ronments, the upper mass cutoff is increasing with environ- 
mental density. The best-fit dividing line from S I3.2| provides 
a reasonable color division between the two sequences in all 
environments as shown by the dash-and-dotted lines. 

As noted in Paper III, the mean positions of the se- 
quences as shown by the grey regions in Fig. [5] do not vary 
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Figure 10. Colour versus concentration index across different environments for stellar masses around 10 10 A4q. The points and solid 
contours represent galaxies in our sample. The dashed ellipses are fixed in each panel outlining the locus of the red- and blue-sequence 
galaxies. The dominant effect is clearly the increase in the fraction of galaxies on the red sequence with increasing density; while the 
locus of the sequences become marginally redder. The dash-and-dotted lines show the best-fit color divider for the mean mass. 



strongly with environment. This is also evident if one plots 
colour versus concentration index as a function of environ- 
ment for different mass ranges. Figure 1101 shows this for 
galaxies with masses around 1O 1O .A40. The clearly domi- 
nant effect is the increase in the fraction of galaxies on the 
red sequence as a function of environment. 



3.5 Environmental and stellar mass dependence of 
the relative red- and blue-sequence numbers 

In lBaldrv et"aT] (l2004bi) . it was shown that the fraction of 
galaxies in the red sequence could be unified by a function 
of environmental density and luminosity. Here, we proceed 
with a similar analysis using stellar mass in place of lumi- 
nosity. The population was divided into 12 environment bins 
and 13 mass bins. The fraction is defined as galaxies redder 
than the best-fit dividing line [Fig. |7{c) and [5]. Figure ITTI 
shows the fraction of red-sequence galaxies (a) versus envi- 
ronment for different stellar masses and (b) versus stellar 
mass for different environments. Qualitatively the plots are 
similar: with increasing red fraction versus stellar mass or 
environment and a greater variation at the lower end of each 
scale. 

Both stellar mass and environment affect the probabil- 



ity of a galaxy being in the red sequence. First of all, we 
found an empirical relation that is based on the sum of pro- 
jected environmental density and stellar mass: 



f r = 0.5 + oi log 10 (E/a 2 + M/az) 



(9) 



(with the constraint ^ f T ^ 1) where f r is the fraction 
of red-sequence galaxies and 01,2,3 are the parameters to 
determine. Secondly, we found an empirical relation more 
naturally related to probability theory that is given by: 



fr 



exp{- 



[(E/6i) b2 + (X/63) 64 ]} 



(10) 



where 61,2,3,4 are the parameters to determine. This function 
has the useful feature of asymptotically tending to and 1 
at the extremes. Figure 1121 shows these empirical relations 
that were determined by minimising \ 2 ■ The values of the 
parameters are (a) 0.438, 10 46 Mpc- 2 , 10 10 - 34 M Q and (b) 
10 a91 Mpc~ 2 , 0.69, 1O 1O 72 X , 0.59. [Note an empirical law 
of the form f r = F(clog E + log M) does not work; the \ 2 
is about 200 higher for a quadratic function F.] 

Both empirical relations provide a reasonable fit to the 
data. The relation in Fig. ll2f b'l provides a marginally better 
fit relative to Fig. ll2f a'l despite the deviation at low masses. 
In addition, the low mass end is uncertain because of the 
difficulty in extracting a reliable low-redshift sample; there 
are systematic uncertainties mostly related to deblending, 
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Figure 11. Fraction of red-sequence galaxies versus environment and versus stellar mass. In panel (a), the symbols and lines represent 
different stellar masses as shown in the legend (log M from 9.0 to 11.6). In panel (b), the lines represent different environmental densities. 
Systematic errors of 0.03 were added in quadrature to the Poisson errors. Note the similarity between the two plots leads to the unification 
schemes shown in Fig. 1121 




(E/10 - 46 Mpc- 2 ) + (M/I0 l0 - M M so]ar ) (E/10 - 91 Mpc- 2 ) - 69 + (M/10 i0 - 7S M. solar )° ™ 

Figure 12. Fraction of red-sequence galaxies versus unified quantities: (a) using Eq. Eland (b) using Eq.EO The symbols and error bars 
represent different stellar masses as shown in Fig. Illf ai. The solid black lines represent the functions shown in the bottom-right corner 
of each panel. 

for example, so me low-redshift system s are in fact deblends better fit. GSMFs for the red- and blue-sequences based on 
of large galaxies dBlanton et al . 2005b). For log A4 > 9.4, the these empirical laws are shown in the Appendix, 
results are more reliable because the average redshifts are 
higher. If \ 2 is determined only using galaxies with log M > 
9.4 then the relation in Fig. 1121 bl provides a significantly 
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4 DISCUSSION 

The advent of large-volume redshift surveys has greatly ex- 
panded the studies of galaxy populations as a function of 
envir onment that ha d traditionally been based around clus- 
ters iDresslerl IT98C0 to the full range from 'void' to clus- 
ter. There are many studies including Paper II to this pa- 
per f § 11.11 that focus on the variation of galaxy proper- 
ties w i th local density <|K auffmann et all2004: |T^naka^^lJ 
| 2004|: iKelm et alJ 20051: lAlonso et alj hood lHaines etldl 
hood: iMateus et al 112009 : IPatiri et allhood . Other studies 
have anal ysed the average density as a function of galaxy 
property jflogg et alJl2003t iBlanton et alj|2005al). the clus- 
terin g properties of galaxies iWild et alJl2005t IZehavi et alJ 
120051: iLi et alJlioOrj . and the relationship between galax- 
ies a nd dark matter derived from galaxy-galaxy weak lens- 
ing jGrav et al] 12004 iMandelbaum et alJhood) . The data 
available for 2 < 0.2 is the most comprehensive, e.g., from 
SDSS and 2dFGRS dat a, but gains have been made at higher 
redshift JWilman et alJl200d lYee et~al]l2005t ICooper et all 
l200dlHbert et alJl2006ft . While local density measurements 
illuminate various environmental trends in the data, galax- 
ies in semi-analytical models are primarily associated with 
dark-matter halos. The main focus of this discussion is com- 
paring models with the data by 'measuring' the analogous 
quantity to projected galaxy density in the models. 

Our results show that the local galaxy population di- 
vides neatly into two types, and that the fraction of each 
type depends exclusively on a combination of stellar mass 
and local environment. Importantly, the environmental de- 
pendence is not a second-order effect, but is at least as im- 
portant as stellar mass in determining the fraction of red 
galaxies in a population. By contrast, in simple galaxy for- 
mation models, the characteristics of a galaxy are usually 
determined primarily by the mass of its dark matter halo, 
which is closely related to the stellar mass of the galaxy. 
That is because the formation time, cooling rate and merg- 
ing history of a halo are most strongly related to its present 
day mass. The most important environmental consideration 
in most of these models is that galaxies at the centre of a 
dark matter halo are tr eated differently from non-central, 
or satellite galaxies (e.g. lCole et al.ll2000l) . Only the central 
galaxies are assumed to possess a halo of hot gas that can po- 
tentially cool and replenish the disk. In addition, there are 
second-order effects related to the large-scale density field 
that, qualitatively at least , can mimic some of the observed 
trends with environment llMaulbetsch et alll200dl . depend- 
ing on how the details of star formation and feedback are 
treated. 

Our results lead naturally to two questions in this con- 
text. One is, how does our measurement of environment, 
E, relate to the dark matter density field Ap/p? And the 
other is whether or not galaxy formation models that are 
successful in other respects are also able to reproduce the 
observed dependence of galaxy colour on environment. We 
will address these questions using the z = output of the 
Virgo Consortium's Millennium Si mulation. The deta ils of 
the dark matter are described in lSpringel et al] (|2005a), and 
we will comp a re wit h the galaxy form a tion m odels of both 
ICroton et all (l200d) and IBower et all <200db These mod - 
els are improved over earlier efforts (e.g. ICole et"ai] l2000l) 
in that, by including a model of feedback from AGN, they 



are able to better match the observed colour distribution as 
a function of galaxy luminosity. Specifically, the inclusion 
of AGN feedback removes most of the bright blue galaxies, 
and increases the col our difference between the red and blue 
populations (see also ISpringel et al]l2005fj) . 

Although both models include feedback from radio-jets 
and AGN ; the models use diffe rent schemes to implement the 
feedback. ICroton et al] ]200d) compute an energy feedback 
rate based on a semi-empirical model involving the mass of 
the host halo and that of the central black hole. Their pa- 
per discusses how the expression they use is related to the 
Bondi accretion rate. In contrast, the IBower et al] ]200d) 
model assumes that AGN feedback will be self-regulating 
if the cooling time is long compared to the sound crossing 
time of the system so that the coo ling of gas tak es place from 
a quasi-hydrostatic atmosphere teinnevil200l) . The model 
also requires that central black hole is sufficiently massive to 
provide heating rate. In addition to the treatment of AGN 
jets, the models also differ in many details of the implemen- 
tation of cooling galaxy merging, star bursts and many other 
factors. 



4.1 Galaxy versus dark matter densities in the 
models 

To construct a mock-observational sample to compare with 
our data, we select galaxies at z = from the simulation. We 
restrict most of the analysis to galaxies with stellar masses 
M > 10 10 Mq, which typically belong to dark matter halos 
with at least 100 particles (so their merger histories are rea- 
sonably well resolved). However, to define the local density 
we will select a luminosity-limited sample (see below) that 
includes lower-mass galaxies. The redshift of each galaxy is 
determined as z = (H Q r z + v z ) /c, where r z and v z are ar- 
bitrarily taken to be the position and peculiar velocity in 
the 2-axis-coordinate of the simulation. The simulation box 
is 714 Mpc on a side, large enough that we can select the 
redshift range 0.01 < z < 0.085 to match our data sample, 
without the need to tile the box. 

We measure E in the simulation using a similar prescrip- 
tion to that used for the data. Specifically, for each galaxy we 
define a DDP as all galaxies in the galaxy catalogue brighter 
than M r = —20, and within cz = 1000 km/s. The distances 
to the fourth- and fifth-nearest neighbours within this DDP 
are calculated for each galaxy in the simulation, and the 
projected density E is computed as for the data. We do not 
attempt to model the SDSS spectroscopic completeness, and 
thus ignore the small effect of only having photometric red- 
shift estimates in incomplete regions (galaxies with large 
uncertainties are not included in the data analysis in any 
case) . 

Figure [T31 shows how the measurement of E compares 



5 Distances, absolute magnitudes and masses in the simulation 
are converted to a cosmology with Ho = 70 km s _1 Mpc — 1 . 
The SDSS r-band luminosity is computed from the B- and V- 
band magnitudes in the catalogue via the transformation r = 
V - 0.42(B - V) + 0.21. This is based on the calibration of 
Ijester et al] 120051) . with an additional 0.1 mag shift to ensure 
the space density of galaxies with M r = — 20 matches that of the 
data. 
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Figure 13. The correlation between dark matter overdensity 
Ap/p, smoothed with a Gaussian kernel of radius 1.78 Mpc, 
and projected gala xy surface density E, from the model of 
ICroton et aTlfcOOrJ) . The data are binned by 0.1 in the logarithm 
of each quantity, and the contours are logarithmically spaced 
by 0.4 with the lowest contour representing 75 galaxies per bin. 
The correlation demonstrates that the observationally defined E 
broadly traces the underlying dark matter density distribution. 

with the smoothed dark matter overdensity, Ap/p, com- 
puted with a Gaussian kern e l of r adius 1.78 Mpc, for all 
galaxies in the ICroton et al.l i200(tl catalogue with M r ^ 
—20. There is a clear correlation, in that galaxies with higher 
S tend to lie in regions that are overdense in dark matter. 
This is encouraging but not surprising. Figure tTH shows how 
this correlation depends on halo mass, by which we mean 
the mass of the largest dark matter halo associated with the 
galaxy (i.e., not the subhalo in which the galaxy may reside). 
The three sets of contours correspond to halos with masses 
within 0.1 of log{M/M Q ) = 14.75, 13.75, 12.75 and 11.75. 
For halos of a given mass, there is a large range in E; thus, 
this estimate of local density is sensitive to variations in en- 
vironment within a given halo. Moreover, the distribution 
of E is only weakly dependent on mass, for M > 10 13 A4©, 
which reflects the fact that the density structure of dark 
matter halos are nearly self-similar. On the other hand, the 
dark matter overdensity on 1.78 Mpc scales is more sensi- 
tive to halo mass, because the size of the smoothing kernel is 
comparable to, or larger than, the typical halo. For low-mass 
halos, where the number of bright M T < — 20 galaxies per 
halo is small, E traces dark matter overdensity reasonably 
well, because both are probing similar scales. 



4.2 Comparison between models and data 

Recently, IWeinmann et al l j2006T) have dem onstrated that 
although the models of Icrotoiiet alJ (I2006T) do a good job 
of reproducing the overall luminosity- dependence of the ob- 
served galaxy colour distribution, they do not succeed in 
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Figure 14. As Fig. 1131 but for four distinct populations divided 
by the mass of the largest dark matter halo in which they are 
embedded. The lowest contour level represents 10 galaxies per 
bin (0.1 in the logarithm of E and Ap/p), and the contours are 
logarithmically spaced by 0.4. The solid contours are restricted to 
galaxies with log Mhalo/MQ = 14.75 ± 0.1; the dotted contours 
represent log Mhalo = 13.75, the dashed contours log.Mh a ; = 
12.75, and the dash-and-dotted contours log M^alo = 11-75. For 
halos of fixed mass, E spans a wide range, as it is sensitive to the 
local environment of galaxies. 



matching the colour distribution in groups and clusters. 
Specifically, they found that the model predicts too many 
faint, red galaxies in these large halos. We will investigate 
this further by characterising environment using the contin- 
uous variable E, which is more closely related to the large- 
scale dark matter density field than to the mass of the halo. 

To divide the m odel into red and blue g alaxies we fol- 
low the approach of IWeinmann et all <l2006f) and make a 
magnitude-dependent colour cut. We find that a separation 
at (B-R) = -0.05 (M r + 20) + 1.83 provides a good separa- 
tion between the red and blue populations of the model (our 
results below are not sensitive to this cut, because the model 
populations are so well separated in colour). Figure [l31 shows 
the fraction of red galaxies, using this criterion, as a func- 
tion of stellar mass and local density E. Qualitatively, the 
models reproduce the observations well, as the red fraction 
increases with both stellar mass and E. However, quantita- 
tively there are some interesting differences. In particular, 
the dependence in the model is actually too strong relative 
to the data. At low densities, the model predicts too few 
red galaxies, with a red fraction that varies from 0.05 to 
0.6 over the plotted mass range. By contrast, over the same 
mass range, the data show a red fraction ranging from 0.4 to 
1.0. Another point of interest is that the mass dependence of 
the red fraction in the models is weak below ~ 10 10,6 A4q, 
while the data show a continuous decrease in red fraction 
with stellar mass over another order of magnitude. At the 
high-density end, the models seem to do a better job of 
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Figure 15. Left panel: The fraction of red galaxies in the ICroton et all 120061) model, as denned in the text, as a function of stellar 
mass and environment. The green, red and black lines represent the model. They are analogous to those shown in Fig. Illf a'l but only 
extend to masses as low as 1O 1O A^0. For reference, the grey lines represent f r from the data parameterised using the unified relation 
of Fig. I12f b - ) for \ogM = 10.0 to 11.6. Right panel: The model galaxies are divided into central galaxies (solid points) and satellite 
galaxies (crosses). 



reproducing the high red fraction and weak stellar-mass de- 
pendence seen in the observations. 

The main source of the environmental dependence in 
the model is simply illustrated in the right panel of Fig. 1151 
where we divide th e sample into cen t ral an d satellite galax- 
ies. As shown by IWeinmann et all J2OO6), most satellite 
galaxies are red, independent of luminosity and density, in 
the models. The central galaxies, on the other hand, have 
a strong stellar mass dependence, and a dependence on en- 
vironment that is considerably weaker than that shown by 
the population as a whole, in the left panel. The overall en- 
vironmental trend is therefore largely due to the fact that 
low-density regions are dominated by central galaxies, while 
the high-density end is dominated by satellite galaxies. The 
main source of the discrepancy with the observations, there- 
fore, can be identified as the low fraction of red galaxies 
amongst the central galaxies in low-density environments. 

Figure 1161 shows the analog of Fig. 1151 for the 



iBower et alJ i2006T) models. In these models, the 
colour- lumino s ity r elation is significantly different from 
ICroton et alJ (120061) . so we have had to adopt a differ- 
ent definition of red galaxy, with a colour cut given by 
(B - R) = -0.07 (M r + 20) + 1.25. For this model, the 
agreement with the observations is remarkably good. The 
difference lies in the colour distribution of the central 
galaxies, for which, at all dens ities, a larger fraction lie 
in the red sequence than in the ICroton et al] model. This 
is not surprising, as the main difference between the two 
models lies in the treatment of AGN feedback, within 
central galaxies. Therefore, the satellite populations are 
largely unaffected. 

4.3 Implications 

IWejnmann et all ((2006) showed that the ICroton et al] 

( 2006) model fails to predict the presence of faint blue galax- 
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Figure 16. As Fig, 1151 but with the galaxy formation model of iBower et alj 1200(1) . The main difference is that the central galaxies 
have a higher red fraction in all environments, leading to better agreement with the observed, overall environmental trend. The model 
parameters were not adjusted to match these data. 



ies in groups and clusters, and link this discrepancy to the 
unifo rmly red sat ellite population that is s een here in both 
the ICroton et alJ and IBower et alJ (|2006) model. This is 
quite a different effect from the discrepancy revealed by our 
results, as our data show that the fraction of red galaxies 
in the densest environments is nearly as high as the models 
predict. Instead, our a nalysis shows a problem with the op- 
posite population in the lCroton et al J model, in that the cen- 
tral galaxies in low-density environments are too frequently 
blue. Interest ingly, this pr oblem does not appear to be as 
severe m 

theH 

ower et all model. However, it is not clear 
whether this improvement is due to their different prescrip- 
tion for AGN feedback, or to some other aspect of the model. 
Unambiguously identifying the root cause of this difference 
requires a more detailed comparison between the two mod- 
els, which is beyond the scope of this paper. 

The comparison between models and data that we 
have discussed and shown in Figs. H5Hl(il is only the red- 
blue sequence galaxy fraction; clearly colour-mass, colour- 
concentration relations and GSMFs could also be compared 



as a function of environment. However, colour-mass and 
colour-concentration relations are more dependent on the 
detailed physics of stellar and galaxy evolution such as 
chemical and dynamical evolution; whereas the red fraction 
(based on a colour cut dividing the bimodal population) rep- 
resents the the most significant variation with environment. 
The focus of this paper is on the unified relation of the red 
fraction versus stellar mass and environment, and we have 
shown that semi-analytical models with their distinction be- 
tween central and satellite galaxies can match this with an 
appropriate feedback d escription (Fi g. [lCf. Of course this 
does not mean that the lBower et al.| (J2006) model has the 
right answer for the galaxy bimodality as there are many 
other pieces of data to match, e.g., lower mass galaxies, 
colour relations and GSMFs. 

Given that the IBower et al. (2006) model produces a 
good match to the unified relation we are now in a position 
to interpret it in the context of the model. In the model, the 
likelihood that a galaxy lies on the red sequence depends 
on two factors. Firstly, the probability that the AGN feed- 
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back is able to prevent the collapse of further gas depends 
on the mass and dynamical time of the halo, and the mass 
of the central black hole. Both of these factors become more 
effective in larger mass galaxies (through the correlation be- 
tween central galaxy mass and halo mass, and through the 
correlation between bulge mass and black hole mass, respec- 
tively) . Secondly, the probability that a galaxy is a satellite 
(and therefore unable to accumulate cold gas) increases as E 
increases. As the galaxy mass increases, the galaxy is more 
likely to be central and red regardless of the E parameter. 
Conversely the galaxies with the lowest mass are very un- 
likely to be central if E is high. Thus the sensitivity to E 
depends strongly on galaxy mass, and this dependence is 
reflected in the unified relation that we have presented. 

The success of the model in representing the impact 
of environment may at first seem surprising since the mod- 
els presented do not yet include many environment-related 
process such as ram pressure stripping of cold disk gas and 
tidal interactions between satellite galaxies. In the models, 
interaction with the environment occurs because the galax- 
ies external gas reservoir is removed as it becomes a satellite 
galaxy orbiting within a larger halo. However, because the 
models assume a rapid interchange between the cold disk 
gas and the external gas reservoir, the loss of the external 
halo quickly suppresses the formation of disk stars. It seems 
that this provides a good approximation to the effect of the 
complete suite of environmental interactions. The caveat, of 
course, is that the process (as currently model) seems overly 
effective so that too few satellite galaxies are able to continue 
star formation, even in if the parent halo has relatively low 
mass. Clearly this is an important area in which the models 
can be developed by introducing a more complete descrip- 
tion for the removal of the external reservoir. 



5 SUMMARY 

We have analysed a primary sample of 151 642 galaxies from 
the SDSS determining stellar masses, rest-frame colours and 
projected environmental densities. The peak, or character- 
istic mass, of the GSMFs increase by about 0.4 dex from 
logE ~ -1.2 to logE ~ +1.3 (Fig. E). The galaxy popu- 
lation is bimodal with a red and blue sequence for which 
the colour-mass and colour-concentration index relations do 
not depend strongly on environment (Figs. 191- 1101 . In con- 
trast, the fraction of galaxies on the red sequence depends 
strongly on both stellar mass and environment (Fig. II IB . 
The stellar-mass and environmental dependence of the red 
fraction can be unified using empirical laws of the form given 
by Eqs. OUCH (Fig. 

In order to compare models with the data, we com- 
puted the equivalent of galaxy projected neighbour density 
for semi-analytical models based on the Millennium Simu- 
lation. This showed that E broadly traced the underlying 
dark- matter density distribution (Figs. I13H141 . The frac- 
tion of red-sequence galaxies versus stel lar mass and envi- 
ronment were com puted for two models IjBower et al . 2006; 
ICroton et al"ll2006h and compared with the unified relation 
('Figs. IT51ll6[l . Both models produ ce qualitative ly the right 
trends with environment; with the lBower et al.l model pro- 
viding a good quantitative match. 
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APPENDIX A: MORE ON THE GALAXY 
STELLAR MASS FUNCTIONS 

Galaxy luminosity functions and GSMFs are typically fitted 
using Schechter functions; including galaxy type-dependent 
functions. There is a logical inconsistency here: the sum 
of Schechter functions do not sum to a Schechter function 
in general. The main results in this paper suggest another 
method of defining type-dependent functions. An overall 
GSMF is defined by a Schechter-based function while the 
red- and blue-sequence functions are defined using a proba- 
bility or fractional function (e.g. Ea. IHlorllOfl. 

Another key point is that an overall luminosity or 
mass function is in many cases not well parameterised by a 
standard Schechter function. Modifications to the Schechter 
funct ion include: adding <j> rea = 10 (a + i ' M) ijMadgwick et alJ 
12003^1 : adding a Gaussian function l|Mercurioet^uT l2006"): 
and using a double Schechter function (Paper I) given by 

^ x AX = e-*'x* [tf (§) Q1 + « (A) Q2 ] |* (Al) 

where <j)x dX is the number density of galaxies with lumi- 
nosity or mass between X and X + dX. One can also specify 
a?2 < ot\ so that the second term will d ominate at the faintest 
magnitudes fcf. iBlanton et al"ll2005b|) . The first term then 
acts to increase the number density around X* . 

The standard Schechter fits are shown in Fig. |5fb) 
over a reduced mass range. In order to fit the mass func- 
tions down to lower masses, we used double Schechter func- 
tions with fixed a?2 = —1.5. The fixed faint-end slope was 
the weighted best- fit average over all e nvironments and 
was consistent with lBlanton et al 1 j2005blV s analysis of low- 
luminosity galaxies. Figure IAT1 shows the fitted GSMFs for 
different environments. The red- and blue-sequence GSMFs 
were then determined by multiplying each overall GSMF by 
/ r (E,A4) and 1 — f r from Eg. 1^1 or 1101 with the parameters 
shown in Fig. The average log E value for each environ- 
ment, the approximate fractional contribution to the stellar 
mass density (/ sm d) and the double Schechter parameters 
(log M* , tf>*, on, 4>2, 012) are shown in each panel. Note with 
the IMF and mass-to-light ratios used in this paper, f^tars ~ 
0.002 (cf.lPersic & Saluccll992HBaldrv fc Glazebrookl2003l: 
iNagamine et al.ll2006h .~ 
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Figure Al. Galaxy stellar mass functions. The dashed lines represent the double Schechter fits to the overall GSMFs. The solid lines 
represent the red- and blue-sequence GSMFs defined by multiplying the overall fit by f r and 1 — f r from Eq. II III while the dotted lines 
represent the sequences using Eq. [5] Each overall GSMF is normalised to a total mass of 10 11 Mq. See §EH1 Fig.|S]and text of Appendix 
for further details. 
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